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Non-Hermitian "PT-symmetric quantum-mechanical Hamiltonians generally ex- 
hibit a phase transition that separates two parametric regions, (i) a region of un- 
broken VT symmetry in which the eigenvalues are all real, and (ii) a region of 
broken VT symmetry in which some of the eigenvalues are complex. This transition 
has recently been observed experimentally in a variety of physical systems. Until 
now, theoretical studies of the VT phase transition have generally been limited to 
one-dimensional models. Here, four nontrivial coupled PT-symmetric Hamiltonians, 
H = \p 2 + \x 2 + \q 2 + \y 2 + igx 2 y, H = \p 2 + \x 2 + \q 2 + y 2 + igx 2 y, H = \p 2 + 
\x 2 + \q 2 + \y 2 + \r 2 + \z 2 + igxyz, and H = \p 2 + \x 2 + \q 2 + y 2 + \r 2 + \z 2 + igxyz 
are examined. Based on extensive numerical studies, this paper conjectures that all 
four models exhibit a phase transition. The transitions are found to occur at g pz 0.1, 
g « 0.04, g « 0.1, and g ~ 0.05. These results suggest that the VT phase transition 
is a robust phenomenon not limited to systems having one degree of freedom. 

PACS numbers: 



I. INTRODUCTION 

We believe that in order to advance the theory of VT quantum mechanics it is crucially 
important to answer the following natural question: Is there a VT phase transition in higher- 
dimensional quantum systems, or is VT quantum mechanics just limited to one-dimensional 
models? 

Recently, there have been some modest attempts to study what happens when two PT- 
symmetric systems are coupled and when a PT-symmetric system is coupled to a con- 
ventionally Hermitian system (see, for example, Ref. [lj), but in these studies only trivial 
matrix and harmonic-oscillator models were considered. In early analytical approaches (see, 
for example, Refs. [21 [3]) some progress was made in calculating the C operator for some 
complicated PT-symmetric Hamiltonians; showing that the C operator exists is equivalent 
to proving that the eigenvalues are real 0]. However, calculating the C operator is difficult, 
and C was only calculated to first order in perturbation theory in Refs. [3]. 

In this paper we report a direct numerical attack on some nontrivial coupled Hamilto- 
nians that were first considered in Ref. [3]. Specifically, we consider the four nontrivial 
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Hamiltonians 

H = \p 2 + \x 2 + \q 2 + \y 2 + z#rr 2 ?/, (1) 

H = \p 2 + \x 2 + |g 2 + ?/ 2 + igx 2 y, (2) 

H = \p 2 + \x 2 + |g 2 + |?/ 2 + \r 2 + ±z 2 + z^t/z, (3) 

H = \v 2 + |^ 2 + |<? 2 + 2/ 2 + |?" 2 + \z 2 + zgxyz, (4) 

where the coupling constants g are real parameters. For these Hamiltonians a convincing 
demonstration that there exists a range of g for which the eigenvalues are all real requires 
the accurate calculation of thousands of eigenvalues. Initial nondefinitive investigations 
using comparatively elementary numerical methods suggested that there might be complex 
eigenvalues for all nonzero values of g [5j. If this were true, then we would be forced to 
admit that VT quantum mechanics is theoretically interesting but of limited scope. 

This paper presents numerical evidence that for each of these complex Hamiltonians there 
are actually ranges of g for which the eigenvalues of these Hamiltonians are all real. Based 
on this numerical work, we conjecture that for H in ([I]) the critical value of g is about 
0.1 and that when \g\ < 0.1 the eigenvalues are all real. For the Hamiltonians in Q, Q, 
and Q the critical values of g are about 0.04, 0.1, and 0.05, and when \g\ is less than 
these values, the eigenvalues are all real. To obtain these results we have used a powerful 
numerical scheme known as the implicitly restarted Arnoldi method [6] . We have calculated 
many thousands of eigenvalues accurate to about one part in 10 6 and have used about 
20, 000 hours of CPU time. Our results suggest that VT quantum mechanics is general and 
robust, and that it extends to genuinely nontrivial higher-dimensional quantum-mechanical 
models. At present we can only draw conclusions based on detailed numerical studies, but 
the results in this paper suggest that VT symmetry may even extend to infinite-dimensional 
quantum-field-theoretic models. 

This paper is organized simply. In Sec. ITT] we describe the VT phase transition and in 



Sec. Ill we summarize our numerical approach and present our results in graphical form. 



Section IV contains brief concluding remarks. 



II. VT PHASE TRANSITION 

A Hamiltonian is VT symmetric if it is invariant under combined space reflection (parity) 
V and time reversal T ■ Such a Hamiltonian is said to have an unbroken VT symmetry if 
its eigenfunctions are also eigenstates of the VT operator. When the VT symmetry of a 
Hamiltonian is unbroken, its eigenvalues are all real even though the Hamiltonian may not 
be Dirac Hermitian [7]. (We use the term Dirac Hermitian to describe a linear operator 
that remains invariant under the combined operations of matrix transposition and complex 
conjugation.) Such a PT-symmetric Hamiltonian is physically relevant because it generates 
unitary time evolution [4J. 

PT-symmetric Hamiltonians usually depend on one or more parameters. The Hamilto- 
nians that have been studied so far typically possess an unbroken-PT-symmetric phase (a 
parametric region of unbroken VT symmetry in which all of the eigenvalues are real) and a 
broken- 'PT-symmetric phase (a parametric region of broken VT symmetry in which some 
of the eigenvalues are complex). The boundary between these two regions is the VT phase 
transition and this transition occurs at a critical value of a parameter. 
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A simple example of a PT-symmetric Hamiltonian that has a VT phase transition is the 
2x2 matrix Hamiltonian [8] 

H ={ r€ s re 1 ")' ( 5 ) 

where the three parameters r, s, and 9 are real. This Hamiltonian is not Dirac Hermitian 
but it is VT invariant, where the parity operator V is 

and the time-reversal operator T is complex conjugation. The region of unbroken VT 
symmetry is s 2 > r 2 sin 2 9. 

Almost all studies of non-Hermitian 'PT-symmetric quantum-mechanical Hamiltonians 
have focused on finite-dimensional matrix Hamiltonians like that in ^ or on one-degree-of- 
freedom Hamiltonians of the type H — p 2 + V(x), for which the condition of VT symmetry 
is V*(x) = V(—x). The Schrodinger eigenvalue problem Hip = Eip for such Hamiltonians 
takes the form of the ordinary differential equation —ip"(x) + V(x)ip(x) — Eip(x). 

The first PT-symmetric Hamiltonian that was studied in detail has the form 

H = p 2 + x 2 {ix) £ , (7) 

where £ is a real parameter [TP] . Dorey, Dunning, and Tateo proved that the eigenvalues 
of the corresponding Schrodinger eigenvalue equation 

- i{j"(x) + x 2 {ix) £ ip(x) = Etp(x) (8) 

are all real when e > [HJ Q2] . These eigenvalues are plotted as functions of e in Fig. [I) The 
critical value of e is 0, and the region of broken VT symmetry is s < 0. As e approaches 
through negative values, complex-conjugate pairs of eigenvalues become degenerate as they 
emerge from the complex plane and then split into real pairs of eigenvalues. The values 
of e at which the eigenvalues become degenerate are sometimes called exceptional points 
[13]. Observe from Fig. [ljthat e = is the limit point of a sequence of exceptional points. 
For this model the decomplexification process is a high-energy phenomenon; that is, as e 
approaches from below, sequentially higher eigenvalues (rather than lower eigenvalues) 
become degenerate. This is because £ is a singular perturbation parameter [Hj. (The 
same kind of limiting process was discovered many years ago for the case of the anharmonic 
oscillator H = p 2 + x 2 + gx 4 , except that for the anharmonic oscillator the exceptional points 
lie in the complex-^ plane rather than on the real-g axis [T5].) 

Many other PT-symmetric Hamiltonians that exhibit VT phase transitions have been 
studied. For example, the Hamiltonian 

H = p 2 + x 4 + iAx, (9) 

where A is a real parameter, exhibits a VT phase transition at A — ±3.169 [16j. In contrast 
to the Hamiltonian as A approaches the phase transition through the region \A\ > 3.169 
of broken VT symmetry, the eigenvalues emerge from the complex domain starting with the 
highest-lying eigenvalues and ending with the lowest-lying eigenvalues. The exceptional 
points remain well separated and do not converge to a limit point because A is a regular 
perturbation parameter |14j . 
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FIG. 1: Real eigenvalues of the "PT-symmetric Hamiltonian plotted as functions of e. The region 
of unbroken VT symmetry is e > 0. In the region e < of broken VT symmetry, the eigenvalues 
emerge from the complex plane at exceptional points as degenerate pairs; these pairs then split into 
real eigenvalues. As e approaches the critical point at from below, this decomplexification process 
begins with the low-lying eigenvalues and terminates at e — at the highest-lying eigenvalues. 



Another example of a PT-symmetric Hamiltonian having a VT phase transition is 

d 2 
cW 



H = +igcos9, (10) 



where g is a real parameter [T7j. Here, the Schrodinger eigenvalue problem is posed on a 
finite rather than on an infinite domain. If the eigenfunctions are required to be 2n periodic 
and odd in 9, the region of unbroken VT symmetry is \g\ < 3.4645 [T7] , 

A physically motivated example of a PT-symmetric Hamiltonian that exhibits a VT 
phase transition was discussed by Rubinstein, Sternberg, and Ma |18j . Their Hamiltonian 
arises in the context of superconducting wires. Again, the Schrodinger eigenvalue problem 

— tJj"(x) — igxTp(x) = Etp(x) (11) 

is posed on the finite domain \x\ < 1 and the boundary conditions are ^(±1) = 0. The region 
of unbroken VT symmetry is g < 12.31, and the behavior of the eigenvalues is qualitatively 



similar to that of H in (10). (Compare Fig. 1 of Ref. [TT] with Fig. 1 of Ref. [IE].) 

The VT phase transition has been seen repeatedly in laboratory experiments. It was first 
observed in optical wave guides [T9ti23] . but it has also been observed in atomic diffusion 
[21] , lasers |25l |26] , superconducting wires [[8] , nuclear magnetic resonance [27] , and most 
recently in microwave cavities [28j and in electronic circuits [29j . The VT phase transition 
is a clear and prominent effect and not a subtle phenomenon; laboratory measurements and 
theoretical predictions have agreed with virtually no error. 

Theoretical studies of eigenvalues (as in Fig. [T]) and experimental observations suggest 
that the VT phase transition is quite generic, but almost all of the theoretical and exper- 
imental work that has been done so far has concerned systems that have just one degree 
of freedom. Indeed, the proof of the reality of the eigenvalues in Refs. [HI [12] relies on 
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establishing a correspondence, known as the ODE/IM correspondence , between ordinary 
differential equations and integrable models. We emphasize that this correspondence in- 
volves ordinary and not partial differential equations. Thus, we are motivated to study in 
Sec. Ill the eigenvalues of the four multidimensional Hamiltonians in (fTl) - (J4|. 



III. NUMERICAL CALCULATION OF EIGENVALUES 

To find the eigenvalues of the Hamiltonians in ([I]) - Q , we follow a calculational approach 
that was used in Ref. [10J. We express the Hamiltonian in harmonic-oscillator-basis states 
[see Eq. (3.2) in Ref. |1U]]. truncate the Hamiltonian matrix to an TV x N numerical array, 
and then calculate the eigenvalues. In this paper the eigenvalues are calculated using the 
implicitly restarted Arnoldi method as implemented by ARPACK [30]. This numerical 
technique works particularly well because the array is sparse] that is, only a small percentage 
of the matrix elements are nonzero. Furthermore, the matrix becomes increasingly sparse 
with increasing dimension. 

As is shown in Fig. 12 of Ref. [TU] . the convergence of the eigenvalues as N increases is 
noisy at first, but starting with the lowest energies, the eigenvalues settle down one- by-one 
to their correct limiting values. In Ref. [10] the eigenvalues of the Hamiltonian H — p 2 + ix 3 
are calculated as the dimension of the matrix ranges up to N — 15. Here, to illustrate the 
convergence of the eigenvalues, we plot in Fig. [2] the eigenvalues of 

H=p 2 + x 2 + ix 3 (12) 

for TV ranging from 2 to 100. Note that the eigenvalues are converging to the values shown 
in Fig. [I] for e = 1. 



A. Calculation of the Eigenvalues of the Hamiltonian ([T]) 

If we use the same approach for H in ([!]) , the numerical convergence is faster than that 
for H in (12) because the matrix is more sparse; only about 5% of the matrix elements are 
nonzero. Using 100 2 x 100 2 matrices (10 8 matrix elements), we calculate the eigenvalues 
for g ranging from to 0.4 in steps of 0.0005. There are 10 4 eigenvalues, but we limit our 
attention to those eigenvalues that have settled down and are changing by less than one 
part in 10 6 as we increase the size of the matrix from 80 2 x 80 2 to 90 2 x 90 2 to 100 2 x 100 2 . 
In Figs. [3] and [4] we plot the real and imaginary parts of those eigenvalues whose real parts 
range from to 16. A key result is given in Fig. [4] it appears that for g greater than about 
0.1 there are many complex eigenvalues, but that there are no complex eigenvalues when 
g < 0.1. This suggests that there is a critical value of g near 0.1. 

From the eigenvalues in Fig. [4] we then select out just the eigenvalues that have nonzero 
imaginary parts. The real parts of these eigenvalues are shown in Fig. [5] Note that the 
real parts of these eigenvalues increase rapidly as g decreases. Thus, if there is a nonzero 
critical value of g at which a VT phase transition occurs, this transition is a high-energy 
phenomenon, as opposed to the VT phase transition of H in Q (see Fig. [I]). 

If we trace a curve through the left-most points in Fig. [5] we can see that this curve rises 
steeply as g decreases, and we believe that this curve becomes infinite at approximately 
g ~ 0.1. Unfortunately, it is not easy to fit a curve numerically through these points 
because they are not very regular. However, when we repeat these numerical calculations 
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FIG. 2: Numerical convergence of the first seven eigenvalues of the Hamiltonian in (12) using the 
implicitly restarted Arnoldi method. At first, the eigenvalues vary chaotically as the dimension N 
of the matrix increases, but eventually, starting with the lowest energies, they settle down to their 
correct numerical values as shown in Fig. [I] at e = 1. 

for the Hamiltonian ([2]), we find that the corresponding points are more regular, and it is 
indeed possible to fit such a curve. 



B. Calculation of the Eigenvalues of the Hamiltonian ([2]) 

For the Hamiltonian ([2]), we construct Figs. [6] and [7J which are the analogs of Figs. [3] and 
[4} Note that the eigenvalues no longer show a breaking of degeneracy as g increases from 0. 
In Fig. [7] we plot the imaginary parts of those eigenvalues. Note that there are no complex 
eigenvalues when g is below about 0.08. 

From the eigenvalues in Fig. [6] we select out just the eigenvalues that have nonzero 
imaginary parts. The real parts of these eigenvalues are shown in Fig. [8| As in the case of 
Fig. |5j these real parts increase rapidly as g decreases. 

Because the left-most points in Fig. [8] are quite regular, we can numerically fit a curve 
that passes approximately through these points. If we seek a curve of the form 

f(g) = a(g - 6) c , (13) 

we find that 

a = 2.32 ±0.18, b = 0.046 ± 0.002, c = -0.615 ± 0.033. (14) 
If we then try a more elaborate curve of the form 

f(g) = a(g-by[-log(g-b)] d , 
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FIG. 3: Real parts of the eigenvalues of H in ([!]) that have converged for a 100 2 x 100 2 matrix 
plotted versus g. The coupling constant g ranges from to 0.4 in steps of 0.0005, and the eigenvalues 
shown have real parts less than 16. The plot is complicated because when g — the eigenvalues 
become increasingly degenerate with increasing energy; there is one eigenvalue of energy 1, two 
eigenvalues of energy 2, three eigenvalues of energy 3, and so on. As g increases, the degeneracy is 
broken. Eventually, one can see energy levels crossing, but most crossings only represent accidental 
degeneracies. One must look carefully to see the exceptional points, which are not so easy to see 
as those in Fig. [Tj In the range of g shown, the lowest-energy exceptional point occurs at about 
g = 0.364, where one of the sixth and one of the seventh eigenvalues become degenerate and form 
a complex-conjugate pair. 

we find that the value of c is consistent with 0, which suggests an improved fitting curve of 
the form 

/(<?) = a[-lo S (g-b)] d . (15) 

The numerical fit then gives 

a = 2.17 ±0.11, b = 0.054 ±0.001, d= 1.67 ±0.5. (16) 

Thus, based on these extrapolations, we are moderately confident that there is a critical 
value of g near 0.05, which is somewhat smaller than the rough estimate of g CT [ t obtained by 
inspection of Fig. |7[ 



C. Calculation of the Eigenvalues of the Hamiltonian ([3]) 

Next we consider the Hamiltonian H in Q . We diagonalize a 30 3 x 30 3 matrix represen- 
tation of H and plot the results in Figs. [9] and [10[ Figure [10| suggests that there is a VT 
phase transition near g ~ 0.25. These eigenvalues have stopped changing as the size of the 
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FIG. 4: Imaginary parts of the eigenvalues whose real parts are shown in Fig. [3j Note that the 
imaginary parts of the eigenvalues are much smaller than the real parts by roughly a factor of 10. 
For the size of the matrix studied, it appears that the critical point is at g ~ 0.1. Below this point 
there are no complex eigenvalues whose real parts are less than 16. 




FIG. 5: Real parts of the eigenvalues in Fig. [4] whose imaginary parts are nonzero. Note that the 
real parts of these eigenvalues grow with decreasing g. We can obtain a more accurate estimate of 
the critical value of g by tracing a curve through the left-most points on the graph. 
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FIG. 6: Real parts of the eigenvalues of the "PT-symmetric Hamiltonian ^ plotted as functions 
of g for < g < 0.4 in steps of 0.0005. All eigenvalues (both real and complex) whose real parts 
are less than 18 are shown. This graph is the analog of Fig. [3j 



matrix size increases from 20 3 x 20 3 to 25 3 x 25 3 to 30 3 x 30 3 . Then, in Fig. 11 we isolate 



just those eigenvalues whose nonzero imaginary parts are shown in Fig. 10 



D. Calculation of the Eigenvalues of the Hamiltonian ([4]) 



The analogs of Figs. [9] and 10 for the Hamiltonian Q are Figs. 12 and 13, However, 
because the eigenvalue degeneracies have been lifted, the left-most points in Fig 



regular enough to perform fits of the form in (13) and (15). The first fit gives the value 



14 



are 



0.049 ±0.001 for the critical point and the second fit gives a value of 0.064 ±0.001 for the 
critical point. Thus, we may estimate that there is a critical value of g near 0.057, which is 



somewhat smaller than what one would guess from Fig. 13 



IV. SUMMARY 



We have presented extensive numerical evidence in this paper that suggests that the 
eigenvalues of the multidimensional trilinear Hamiltonians ([l}j4]) are entirely real when the 
coupling constant g is less than a critical value. If this numerical evidence stands up to 
further study, then we can conclude that these quantum theories possess an unbroken PT- 
symmetric phase for sufficiently small g and display a transition to a broken phase as g 
increases. This would suggest that multidimensional PT-symmetric quantum systems can 
exhibit the same kind of phase transition that one-dimensional quantum systems are known 
to exhibit. On the basis of this work we are tempted to conjecture that even PT-symmetric 
quantum field theories such as an igcj) 3 theory might exhibit a VT phase transition as the 
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FIG. 7: Imaginary parts of the eigenvalues whose real parts are shown in Fig. [6j Note that the 
imaginary parts of the eigenvalues are smaller than the real parts by a factor of roughly 10. For the 
size of the matrix studied, it appears that the critical point is near g « 0.08. Below this point we 
see no complex eigenvalues whose real parts are less than 18. Further analysis (see Fig. [8]) suggests 
that the critical point is near 0.05. 

coupling constant g increases from 0. 

In this paper we have only considered cubic and trilinear interactions because the 
harmonic-oscillator-basis functions that we have used to construct the large numerical ma- 
trices whose eigenvalues we have obtained numerically have exponentially vanishing asymp- 
totic behaviors on the real axis. It would be interesting to study quartic Hamiltonians using 
similar techniques, but we anticipate that the numerical analysis would be much more diffi- 
cult because the harmonic-oscillator-basis functions and the exact eigenfunctions have very 
different asymptotic behaviors on the real axis. 

Finally, we mention again that if there really is a PT-symmetric phase transition for 
the models that we have studied, this transition is generically a high-energy phenomenon. 
Unfortunately, this presents a severe problem for the numerical techniques used in this paper 
because the Arnoldi method can only give detailed information about converged eigenvalues 
for low energies as the size of the matrix is increased. Thus, the principal finding in this 
paper, namely, that the models we have studied exhibit a VT phase transition, can only be 
conjectural in nature. 
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FIG. 8: Real parts of the eigenvalues in Fig. [6] whose imaginary parts are nonzero; the real parts 
of these eigenvalues grow with decreasing g. We can obtain a more accurate determination of the 
critical value of g by fitting a curve through the left-most points on the graph, which are much 
more regular than the left-most points in the analogous Fig. [5] 

of its resources. 
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FIG. 11: Real parts of the eigenvalues in Fig. [9] whose imaginary parts are nonzero. Note that 
the real parts of these eigenvalues grow with decreasing g. We can try to obtain a more accurate 
determination of the critical value of g by tracing a curve through the left-most points on the graph 
and estimating that this curve blows up near g « 0.2. Unfortunately these points are not regular 
enough to do a numerical fit to such a curve. 




FIG. 12: Real parts of the eigenvalues of the 'PT-symmetric Hamiltonian Q plotted as functions 
of g for g ranging from to 0.4 in steps of 0.0005. All eigenvalues (both real and complex) whose 
real parts are less than 18 are shown. 
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FIG. 13: Imaginary parts of the eigenvalues whose real parts are shown in Fig. 12 There appears 
to be a critical point near g « 0.09. Below this point there are no complex eigenvalues in the 



range studied. Further analysis (see Fig. 14) suggests that the critical value of g is near 0.057. 



The eigenvalue behavior shown in this figure is similar to that shown in Figs. [4| [7| and 10 for the 
Hamiltonians in (fTT-[3|. 




FIG. 14: Real parts of the eigenvalues in Fig. 12 whose imaginary parts are nonzero. These 
eigenvalues grow with decreasing g. A numerical fit to the left-most points on the graph predicts 
that the critical value of g is about 0.05. 



